arXiv:astro-ph/0503135vl 7 Mar 2005 


Submitted to ApJ. 

Preprint typeset using I^-T^rjX style emulateapj v. 6/22/04 


SIMULATIONS OF DUST IN INTERACTING GALAXIES I: DUST ATTENUATION 

Patrik Jonsson 

Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064 


T. J. Cox 1 and Joel R. Primack 
Department of Physics, University of California, Santa Cruz, CA 95064 

AND 

Rachel S. Somerville 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore MD 21218 

Submitted to ApJ. 

ABSTRACT 

A new Monte-Carlo radiative-transfer code, Sunrise , is used in conjunction with hydrodynamic 
simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations 
are in good agreement with observations of dust absorption in starburst galaxies, and the dust has 
a profound effect on their appearance. The dust attenuation increases with luminosity such that at 
peak luminosities ~ 90% of the bolometric luminosity is absorbed by dust. In general, the detailed 
appearance of the merging event depends on the stage of the merger and the geometry of the encounter. 
The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be 
predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and 
metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and 
metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust 
luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be 
predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust 
attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier 
studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy 
formation. 

Subject headings: dust — radiative transfer — galaxies: interactions — galaxies: starburst — methods: 
numerical 


1. INTRODUCTION 

Galaxy mergers produce some of the most spectacu¬ 
lar events in the universe. Locally, they are responsible 
for the most lumi nous objects, Ultralumino us Infrared 
Galaxies lULIRGs. [Sanders fc Mirabelll99fi l. At higher 
redsliifts they may be responsi ble for the sources seen 
in the submillimeter llSmall et a, bill 99 7). and they may 
even be a dominant mode of star formation in the early 
universe j|g 2 fB£Lgng_gtal .1 1200 It lElbaz fc Cesarskvl200.‘ii 
but see lWolf et al.l20041 . Because the most luminous ob¬ 
jects are also generally the most dust-obscured, it was not 
until the launch of IRAS that the existence of ULIRGs 
was discovered. What in the optical appeared to be fairly 
unimpressive, albeit peculiar, galaxies turned out to be 
the brightest infrared sources in the local universe. It is 
now clear that starbursts and dust generally go hand-in- 
hand. The large amounts of gas necessary to fuel a major 
starburst bring with them large column densities of dust, 
obscuring the starburst and reradiating the energy in the 
far infrared. In addition, rapid star formation quickly en¬ 
riches the region with metals, increasing the amount of 
dust further. Including the effects of dust is thus crucial 
when studying these systems. 
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Theoretical studies of dust in galaxies have generally 
not leveraged the power of hydrodynam ic simulations 
but rather used fairly simple geometries. iGordon et alJ 
lt!99l . using a Monte-Carlo radiative-transfer model, 
concluded that Small-Magellanic-Cloud-type (as op¬ 
posed to Milky-Way-type) dust was necessary to fit ob¬ 
servations of starburst galaxies, and that their spectral 
energy distribution (SED) was best replicated by a model 
w here a c lumpy shell of dust surrounds the starburst. 
iFerrara et alJ ( 1999 1 calculated the effect of dust in vari¬ 
ous analytical di sk plus bulg e mode ls, also using a Monte- 
Carlo model. I Silva et all (I1998D and iCharlot fc Fall 
(120001 1 modeled dust effects using a simple two-phase 


model, where young stars are located in dense molecular 
clouds while older stars are extinguished only by the dif¬ 
fuse ISM, and showed that their models could fit spectral 
energy distributions of various observed galaxies. How¬ 
ever, all these models had numerous free parameters and 
gave little guidance on how to choose them, thus limiting 
their predictive capability. 

Notable previous efforts at simulating merging galaxies 
using hydrodyna mic N-bodv codes w it h sta r form ation 
and feed back are iMihos fc Hernauistl ( 1994 ll996 f and 
iSoringel ( 200 01. but these efforts did not consider realis¬ 
tic obser vations of their s i mulations inclu ding the effects 
of dust. iBekki fc Shioval (I2001L l2000albfl did use a sim¬ 
ple model for dust attenuation and reradiation along with 
N-body simulations to investigate major mergers. How¬ 
ever, their modeling, using the “sticky particle” method 
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(lSchwarzlll981ll and with no supernova feedback, did not 
allow them to capture essential features of the hydrody¬ 
namics of merging galaxies. This work builds on these 
earlier works, using a full radiative-transfer model to 
study the effects of dust in the most comprehensive and 
sophisticated suite of major-merger simulations d one so 
far. The radiative-transfer model is described in [fonssonl 
I2001 2 . and det ails about t he hy drodynamic simul ations 
can be found in ICoxI (I2004D 3 and I Cox et alJ (I2005H . 

In this paper, we emphasize one particular aspect of 
these simulations: the presence of a tight correlation be¬ 
tween the fraction of energy absorbed and quantities such 
as the luminosity, metallicity, and mass of the simulated 
systems. These relations should be useful for theoreti¬ 
cal modeling of dust absorption in galaxies, for example 
in semi-analytic models of galaxy formation. This pa¬ 
per is the first in a series describing the effects of dust 
in the simulations, the appearance of the simulations 
and comparisons to observed starburst galaxies (Jonsson 
et al. and Lotz et al., in preparation) includin g using 
new nonparametric measures of morphology llLotz et al.l 
12001 . Other papers based on these simulations concern 
the creation of a ha lo of hot, sh ock-heated gas around 
the merger remnant (ICox et aU2004l) . the star-formation 
properties of minor as well as major merger simulations 
(Cox et al., in preparation), 

2. MODEL DESCRIPTION 

The simulations, performed with the “entropy- 
conserving” version of the GADGET SPH code 
ijSnringel et al.ll200lt iSnringel fc Hernauistl f20021 , con¬ 
sist of encounters of identical copies of the galaxies listed 
in Tabled Because of the very large parameter space, 
the number of simulations needed to exhaustively sam¬ 
ple all possible initial conditions would be prohibitively 
large (on the order of thousands). Instead, the strategy 
used was to define a few reasonably realistic cases based 
on information from observations and cosmological sim¬ 
ulations. 

The main differences from previous work are the in¬ 
clusion of efficient supernova feedback, a correct nor¬ 
malization of the star-formati on law, and an accurate 
treatment of hydrodynamics ( ICox et aD '21)l) .~>i. These 
served to make the merger-induced starbursts less in¬ 
tense. Peak star-formation rates for major mergers of 
typical late-type spiral galaxies are generally in the range 
30 — 50 Mq yr _1 , except for very short spikes. 

2.1. Model Galaxies 

Several different dynamically self-consistent disk- 
galaxy models, listed in Table El are used in this study. 

The Sbc galaxy was modeled using median pr oper- 
ties for Sbc galaxies from [Roberts fc Havnesl 1 1994 1 and 
bulge information from Ide^onel <|199 (t 1L Lower- and 
higher-mass variants of the Sbc galaxy, called Sbc- and 
Sbc+, respectively, were modeled using the 25th and 
75th percentiles of the distributions of properties. An¬ 
other set of galaxies, the “G-series”, was based on mean 
relations from the Sloan Di gital Sky Survey llShen et alJ 
B . gas proper ties from [ Bell et al ] 12003) . and bulge 
information from Ide .Tonel lll99fl These galaxies have 

2 Available at http://sunrise.familjenjonsson.org/thesis 

3 Available at http://physics.ucsc.edu/~tj/work/thesis 


TABLE 1 

The galaxy models used for the merger simulations. 


Model 

M vir a 

(M 0 ) 

M b h 

(M 0 ) 

h c 

trot 4 

(kms -1 ) 

® O 
N N, 

Sbc+ 

9.28 • 

10 11 

1.56 ■ 

10 11 

0.58 

210 

1.12 

Sbc 

8.12 • 

10 11 

1.03 ■ 

10 11 

0.58 

195 

1.00 

G3 

1.16 • 

10 12 

6.22 ■ 

10 10 

0.20 

192 

1.00 

Sbc- 

3.60 • 

10 11 

4.98 ■ 

10 10 

0.58 

155 

0.70 

G2 

5.10 ■ 

10 11 

1.98 ■ 

10 10 

0.24 

139 

0.56 

G1 

2.00 • 

10 11 

7.00 ■ 

10 9 

0.29 

103 

0.40 

GO 

5.10 • 

10 10 

1.60 ■ 

10 9 

0.38 

67 

0.28 


Note. — For the Sbc galaxies, 9 mergers with different 
orbital geometries were used. For the other galaxies, only a 
prograde-prograde encounter was used. 

1 Virial mass. 

2 Baryonic mass. 

3 Gas fraction (of baryonic mass). 

4 Circular velocity. 

5 Metallicity (gas and stellar). 

significantly less gas than the Sbc models and cover a 
much larger range in mass. Detailed information a bout 
these and other model galaxies can be found in Ic^l 
(120041) . The galaxies we r e assi gned a metallicity consis¬ 
tent with iZa.ritskv et akl l|1 9941) . used for both the stellar 
SEDs and the initial gas metallicity. All in all, the seven 
galaxy models cover a range of about 100 in baryonic 
mass, 3 in gas fraction, and 4 in metallicity. 

2.2. Galaxy Merger Setup 

Identical pairs of each of the seven galaxies in Tabled 
were put on a parabolic orbit with the disks prograde. 
One of the disks was in the plane of the orbit, while the 
other was tilted by 30°. Eight additional mergers of the 
Sbc galaxies were also simulated, exploring variations in 
galaxy orientation and encounter orbi t. De tails about 
the merger setup can be found in ICoxI ( 2004 ). 

2.3. Radiative-Transfer Model 

The Monte-Carlo radi ative-transfer code Sunrise , simi¬ 
lar to the DIRTY co d e dCordon et al.[l2001fl . is described 
in detail in lJonssonl (l2004fl and in a future paper. Briefly, 
the geometry of gas and stars along with the detailed 
star-formation histories of the galaxies at roughly 50 dif¬ 
ferent points in time were taken from the N-body simula¬ 
tions and used as inputs to the radiative-transfer calcula¬ 
tion. Outputs are multi-wavelength images of the system 
from a number of viewpoints, as well as luminosity ab¬ 
sorbed by dust and reradiated in the infrared. Sunrise 
is available to interested prospective users. 4 

The stellar SEDs were taken from the Starburst99 
(V4.0. Ibeitherer et al.Nl 99*91) stellar population s ynthesis 
mode l. The IMF used, similar to a Kroupa IMF IlKrounal 
SH, has a high-mass slope of -2.35 from the upper limit 
of 150 M 0 down to 1 Mq. Below this mass it is flat down 
to the lower limit of 0.1 Mq. The disk stars existing at 
the start of the simulation were assumed to have been 
forming at a uniform rate for the previous 8 Gyr. The 
bulge stars were assumed to have formed in an instanta¬ 
neous burst 8 Gyr earlier. 

4 The Sunrise web site is http: //sunrise. f amiljenj onsson.org 
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Time (Gyr) 


Fig. 1. — Overview of a major-merger event. The graph shows 
star-formation rate as a function of time for the prograde-prograde 
merger between two Sbc galaxies. The figures show the distribution 
of gas in the simulations at the indicated points in time from a 
viewpoint perpendicular to the orbital plane. 


To be able to describe the simulation geometry with 
sufficient accuracy while keeping computational require¬ 
ments reasonable, Sunrise uses an adaptive grid. The 
grid refinement criteri a are ad j ustabl e; the parameters 
used are described in iJonsso 3 HU. The number of 
cells in the resulting grids ranged from around 70,000 at 
the start of the simulation to almost 300,000 at the time 
of peak star-formation rate and highest densities. 

The dust grai n model used was the " 11 = 3.1 ” Milky- 
Way model by IWeingartner fc Drain^ JgOOlji with a 
dust-to-metal ratio md/rn m = 0.4 (IPwe kl ll998l) . The 
radiative-transfer calculation was performed for 22 differ¬ 
ent wavelengths between 21nm and 5/im, including the 
Ha and H/3 nebular emission lines. One million Monte- 
Carlo rays were traced for each wavelength. Rather than 
self-consistently computing the infrared dust-em ission 
spectrum, the templates of lDevriendt et all (I1999T) . in¬ 
terpolated to the correct luminosity, were used. 

3. RESULTS 

First, some general results from the radiative-transfer 
simulations will be presented. Then we will continue with 
the specific results concerning dust attenuation in the 
simulations. Throughout this paper, “attenuation” will 
be understood to refer to the fraction by which lumi¬ 
nosity is decreased due to dust, either averaged over all 
directions or along a specific line of sight. The term at¬ 
tenuation is preferred over absorption to emphasize that 
both absorption and scattering processes in a complex 
geometry contribute to a net decrease (or, in rare cases, 
increase) in the emerging radiation. 

3.1. Overview of a Major Merger 



time/Gyr 

Fig. 2. — Time evolution of bolometric luminosity and dust 
attenuation for the prograde-prograde Sbc merger simulation. The 
bolometric luminosity of the system is shown with a dashed line, 
while the amount of luminosity absorbed by dust and reradiated 
in the infrared is shown with a dotted line. The solid lines show 
the luminosity not absorbed by dust, emerging in the UV/visual 
(more precisely the luminosity inferred by an observer under an 
assumption of isotropy), one line for each of the 11 viewing angles. 
The luminosity emerging in the UV/visual is practically indepen¬ 
dent of the bolometric luminosity and shows significant variation 
with viewing angle. 


While the detailed behavior of a galaxy merger de¬ 
pends on the initial setup of the two galaxies, the gen¬ 
eral sequence of events remains the same. Figure^ shows 
the star-formation rate of a merger between two Sbc 
galaxies, along with images of the system at different 
points in time. As the two galaxies initially approach 
each other, they are largely unaffected by each other’s 
presence. Tidal forces disrupt the galaxy disks at the 
time of the first close passage. The star-formation rate 
then increases due to gas being driven to the centers as 
the galaxies separate. The star-formation rate remains 
elevated as the galaxies turn around and again approach. 
There are several short bursts of star formation as the 
two components merge. After coalescence, the merger 
remnant relaxes, and the star-formation rate trails off 
over a timescale of about 1 Gyr. 

The evolution of the bolometric luminosity of the same 
simulation is shown in Figure [2j The bolometric lu¬ 
minosity of the system increases by almost an order of 
magnitude during the course of the merger, closely fol¬ 
lowing the star-formation rate, as expected for a star- 
burst system where young stars dominate the energy 
output. However, the fraction of luminosity which is 
absorbed by dust scales with luminosity in such a way 
that the luminosity not absorbed by dust typically stays 
almost constant. As the star-formation rate and hence 
the luminosity increase, the column densities of gas and 
dust also go up, and with them the attenuation; the 
two factors appear to conspire to keep the UV/visual 
luminosity roughly constant. This effect can be seen 
in Figure |3 where the solid lines, showing the lumi¬ 
nosity escaping in different directions, show little vari- 
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ation with time . This tendency was also observed by 
iBekki fc Shioval l20Q0afl i n their sim ulati ons, a nd in ob¬ 
servational studies ifAdelberger fe Steidell2000l) . 

While there is little temporal variation in the lumi¬ 
nosity not absorbed by dust, there is, in general, a sub¬ 
stantial difference between the flux escaping in different 
directions. Because the initial systems consist of disk 
galaxies, radiation escapes preferentially out of the plane 
of the galaxies. This creates a difference of around a fac¬ 
tor of 2 in the escaping flux, clearly visible in Figure [3 
The merger remnant is more spheroidal and thus the dif¬ 
ference with viewing angle is smaller at the end stage 
of the mergers. As will be shown below, this variation 
in viewing angle means that quantities not dependent 
on direction, such as the fraction of bolometric luminos¬ 
ity which is absorbed, can be predicted with significantly 
smaller scatter than directional quantities like the emerg¬ 
ing luminosity in a specific direction. 

3.2. Simulated Images and Spectra 

As mentioned in the model description, the out¬ 
puts from the radiative-transfer calculations are multi¬ 
wavelength images from a number of viewpoints or, more 
precisely, an SED for each pixel in the images. To fa¬ 
cilitate comparisons with observations, it has been in¬ 
tegrated into images in a number of broadband filters 
(listed in Tabled covering wavelengths from the GALEX 
FUV band to the IRAC2 band on the Spitzer Space Tele¬ 
scope and also collapsed into a spatially integrated SED. 
In total, about 25 major-merger simulations have been 
run through the radiative-transfer code, each at roughly 
50 points in time. For each point in time, images and 
spectra have been generated from 11 viewpoints, equally 
distributed in solid angle, and in 12 different filters. This 
results in a grand total of roughly 10 5 images and 10 4 
spectra. The analysis of these is the subject of a com¬ 
panion paper (Jonsson et al., in preparation), but to il¬ 
lustrate the outputs generated some images are shown in 
Figure|3and a corresponding spectrum in Figure 21° 

3.3. Comparing to Observations 

In order to determine how well our simulations mimic 
actual starburst galaxies, they must be compared to ob- 
servat ions. While this comparison is described in ilonssonl 
lt200l and is the subject of a companion paper, example 
results are included here to illustrate that the simulations 
appear to replicate the properties of observed galaxies. 
After accounting for dust effects, the simulations have 
an absolute r-band magnitude in the range -21.5 to - 
22.5 and a u — r color of 1.3 to 2.2, falling in the region 
of b right, blue galaxies in the Sloan Digital Sky Survey 
llBaldrv et alJl2004 b Even 1 Gyr after the merger, there 
is enough ongoing star formation in the merger remnant 
that it is among the blu e galaxies. In recent simulations, 
ISoringel et all (I2005albfl have shown that feedback from 
an active galactic nucleus can help truncate star forma¬ 
tion, which would let the merger remnant redden more 
quickly. In any case, the inclusion of dust is crucial for 
the agreement; without it, the systems are far too bright 
and blue to agree with observed galaxies. 

5 Movies and the full set of color images of the 
merger simulations can be found on the Internet at 
http://sunrise.familjenjonsson.org/thesis 


The simulations w ere also compared to observations by 
IMeurer et all 1119991 hereafter MHC) and IHeckman et all 
1199811 . Both of these studies looked at correlations be¬ 
tween dust attenuation and other properties of starburst 
galaxies. MHC looked at the relation between dust ab¬ 
sorption, indicated by the far-infrared over ultraviolet 
flux ratio, and the ultraviolet spectral slope. In a sam¬ 
ple of moderately luminous starburst galaxies observed 
by International Ultraviolet Explorer (IUE) and IRAS, 
they found a fairly ti ght correlation b etwee n the two pa¬ 
rameters. However, iGoldader et H3 ll200l examined a 
small sample of more luminous LIRGs and ULIRGs, and 
found that they depart from the relation seen by MHC 
in the sense that their UV color is too blue for their 
infrared luminosity. When the same quantities are ex¬ 
tracted from the simulations, the agreement, shown in 
Figure|Sl is remarkable. The simulations follow the MHC 
relation when similar-luminosity systems are selected. In 
contrast, when the highest-luminosity subsample of the 
simulated galaxies is selected, th ey are found in the sam e 
region as the LIRGs/ULIRGs of lGoldader et al.l (I2002H . 
It should be emphasized that this is not a result of fit¬ 
ting the simulations to these observations, but rather a 
prediction of our initial conditions which were selected, 
a priori, to be re alistic fo r local spir al galaxies. Unlike 
in earlier studies lIGordon et alJlll997 ‘). this agreement is 
contingent on the use of Milky-Way dust. SMC-type dust 
leads to far too red ultraviolet slopes and is inconsistent 
w ith the MHC re l ation i n our simulations. 

iHeckman et all (1998) studied correlations between 
dust absorption, infrared and ultraviolet luminosity, 
metallicity, absolute magnitude, etc. The agreement 
with those results, though not shown here, is also en¬ 
couraging. These, while not exhaustive tests, indicate 
that our simulations provide a good replication of local 
starbursts. A future paper (Jonsson et al., in prepara¬ 
tion) will explore the comparison between simulations 
and observations in detail. 

The simulations can also be used to study the per¬ 
formance of star-formation indicators in the presence of 
dust. While the far-infrared luminosity provides a robust 
indication of the star-formation rate for starburst sys¬ 
tems, both the Ha and the ultraviolet luminosity suffer 
from severe dust effects llJonsso 3Hool. Publis hed dust 
corre ctions based on the Balmer line ratiojCalzett^et^L 
1994) or the ultraviolet spectral slope (iBell fc Kennicutt 
2001) improve the estimate, but the star-formation rate 
can still be underestimated by up to an order of magni¬ 
tude at the highest luminosities, where the dust correc¬ 
tions perform the worst. A detailed analysis of the star- 
formation indicators in the presence of dust is planned 
for a future paper. 

3.4. Dust Attenuation 

Looking at the nine Sbc simulations shown in Figured! 
there is a tight correlation between the bolometric lumi¬ 
nosity of the system and the bolometric dust attenuation, 
averaged over all directions. For a given bolometric lu¬ 
minosity a uniquely determined fraction of luminosity is 
absorbed by the dust; it does not seem to matter if the 
system consists of two barely-interacting disks, a merger- 
driven starburst (on various orbits), or a post-starburst 
remnant. Simulations using galaxies of different mass or 
metallicity follow similar correlations offset to higher or 
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Fig. 3. — Images of the Sbc-Sbc prograde-prograde major merger at four different stages. On the left, about 0.6 Gyr into the simulation, 
the galaxies are passing each other for the first time. A line of induced star formation can be seen in the region between the two galaxies, 
where the gas in the two disks is being shocked. The next image, at 1.6 Gyr, shows the galaxies coming together for the second time. This 
collision is almost perfectly radial, concentrating a lot of gas and dust in the center, as evidenced by the clear dust lane. The third image, 
at 1.7 Gyr, shows the final coalescence of the two galaxies, with characteristic tidal tails. Finally, the rightmost image shows the merger 
remnant 2.0 Gyr after the start of the simulation. While the stellar distribution is largely spheroidal, a large fraction o f the gas not consumed 
in th e starburst has now cooled and reformed a disk, visible as a dust lane. (Disk reformation was also reported bv ISpringel &; HernauisH 
120051 ) The images co ver 67kpc squared a nd are in the SDSS g band. (The electronic edition shows urz color composite images generated 
using the algorithm of ILupton et a DIESS) 



Fig. 4.— The spectral energy distribution associated with the 
second image of the Sbc-Sbc prograde-prograde major merger in 
Figure 0 1.6 Gyr into the simulation. The solid (blue in the elec¬ 
tronic edition) line shows the intrinsic stellar spectrum, while the 
dashed (red) line shows the emerging spectrum after taking dust at¬ 
tenuation into account. At wavelengths longer than 5um. the dust 
e missi on spectrum is taken from the templates of|Devriendt_etjahj 
l. i 99 If) . The two emission lines in the optical are Ha and H/3. In the 
ultraviolet, the dust attenuation is over an order of magnitude, and 
the signature of the well-known “2200A bump” in the Milky-Way 
extinction curve is easily discernable. 


lower attenuation. 

Theoretically, it is expected that luminosity (for star¬ 
bursting systems largely determined by star-formation 
rate) and dust absorption should correlate in the sim¬ 
ulations, as both of these quantities are driven by gas 
density. Concentrating the gas to larger densities will 
increase the star-formation rate, and hence the bolo- 
rnetric luminosity, through the Schmidt law used. Con¬ 
centrating the gas also increases the gas, and hence 
dust, column density, increasing dust absorption. Sev¬ 
eral observational studies have also concluded that dust 
attenuation generally seems to increase with galaxy 


luminosity and star-format i on ra te tWang_^JHeckman| 
1996j_|Adelbercer&Steide|_l2000t_|HogWn^^^^^^^P 

Sullivan et all 120011 iViih et a.l.ll2()03t iBuat et all 12004). 

In the following section, a toy model for estimating dust 
absorption will be presented. This will, in turn, motivate 
the fitting formula used to estimate the dust absorption 
in the simulations. 

3.5. A Toy Model for Dust Attenuation 

In order to come up with a simple model for how 
dust absorption should depend on luminosity, mass, and 
metallicity of a galaxy, consider a constant-density sphere 
of star-forming gas. For a sphere, the density is given by 

p ex M g R~ 3 , (1) 

where M g is the gas mass and R is the radius. According 
to the Schmidt law used to estimate the star-formation 
rate in the simulations, the star-formation rate density 
is 

p* oc p 3/2 . (2) 

The total star-formation rate M* is thus 

M* oc M 3 J 2 R- 3 ' 2 . (3) 

The optical depth of dust in the sphere will depend on 
the column density and the metallicity Z of the gas, 

t oc ZRp oc ZM g R~ 2 . (4) 

Eliminating R using equation [21 we get 

t oc ZM+ /3 M~ 1 oc ZL^Mg 1 , (5) 

where the last proportionality comes from assuming 
that the bolometric luminosity L is proportional to the 
star-formation rate. Finally, once the optical depth is 
determined, the absorbed fraction of luminosity (i.e., 
the attenuation) in a medium where luminous and ab¬ 
sorbing material is uniformly mixed is given by (e.g. 
iCalzetti et alill994 ) 

^ = l-I(l-e- T ). (6) 

L t v ' 

Actually, equation [S] is appropriate for a plane-parallel 
slab, not a sphere. However, the purpose of the toy 
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Fig. 5. — The relation b etween the IR/UV flux ratio and the UV spectr al slope for the Sbcrnerger simulations (shaded region), 
compared to the results from [CTeurer et alJ <119991 hereafter MHC, crosses) and IGoldader et alJ <20021 diamonds/triangles). On the left, 
only simulated galaxies with bolometric luminosity Z/bol < 2 • 10 11 L© have been included. This low-luminosity sample agrees fairly well 
with the MHC correlation, which is for galaxies in this luminosity range. On the right, only the highest-luminosity simulated galaxies, 
with Lbol > 7 • 10 11 L <0 have been includ ed. These points depart completely from the MHC galaxies and instead occupy the region of 
LIRGs/ULIRGs from the lGoldader et alJ (1200211 sample. This agreement was not a result of fitting the models, but rather a prediction 
from our simulations, with no adjustment of parameters. This and other results indicate that our simulations provide a reasonably good 
replication of the properties of local starbursts. 



Fig. 6 .— The relation between bolometric luminosity and at¬ 
tenuation (fraction of bolometric luminosity absorbed by dust, av¬ 
eraged over all directions) for the nine different Sbc-Sbc major 
merger simulations on various orbits. Each point is a simulation 
snapshot. For a given luminosity, the attenuations lie in a narrow 
range, regardless of encounter geometry, stage of merger and the 
fact that different simulations cover different luminosity intervals. 
All simulations were started with solar metallicity gas. (In the elec¬ 
tronic edition, the diff erent c olors identify the different simulations 
in Appendix A of lCoJ2004 'l 


model is to find a simple, physically motivated, fitting 
formula. As will be shown in the next section, equation E] 
describes the behavior of the simulations well. Hence, in 
the interest of simplicity, it will be used. 


This toy model has obvious limitations: It neglects 
scattering, and the exponents 4/3 and —1 depend on 
the assumed geometry. Furthermore, the model is really 
more representative of an individual star-forming region 
than an entire galaxy, so included in the constant of pro¬ 
portionality in equation 0 is the number of such regions 
in the galaxy. If this number depends on the properties 
of the galaxy, a reasonable assumption, it will change 
the dependence on the different quantities in equation 0 
Finally, the attenuation of the bolometric luminosity is 
an average of the attenuation at all wavelengths, and be¬ 
cause the system in general is optically thick at short 
wavelengths and more or less optically thin at longer 
wavelengths, the behavior is more complicated than the 
simple equation El Real galaxies are thus more compli¬ 
cated than our assumptions, but our simple model gives 
a general description of the trends for the effects of dust. 


3.6. Fitting Functions 

While L oc M* is a good approximation for starburst¬ 
ing galaxies, where young stellar populations dominate 
the luminosity, this assumption is not good in general. 
For this reason, simultaneous dependence on both L and 
M* will be retained. 

Guided by the toy model, one would expect the gas 
mass to be the dominant factor in the relation. How¬ 
ever, fits using total baryonic (gas plus stars) mass, Mb, 
instead of gas mass had lower residuals. For this reason, 
baryonic mass is used for the fits. 

Motivated by the above toy model and these consid¬ 
erations, the general fitting function to be used for the 
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TABLE 2 

Fits to the bolometric attenuation. 


TO" 

a a 

/3 a 

7 a 

<5 a 

cr b 

Comment 

1.25 

1.02 

0.41 

0.39 

-0.82 

0.04 

Full fit 

0.93 

1.10 


0.61 

-0.68 

0.04 

Indep. of L 

2.10 

0.91 

1.10 


-1.03 

0.05 

Indep. of M* 

2.09 

0.94 




0.14 

Dep. on Z only. 

2.34 


0.28 



0.12 

Dep. on L only. 

1.70 



0.26 


0.10 

Dep. on M* only. 

2.21 




0.16 

0.16 

Dep. on Mb only. 

0.32 

0.30 

-0.15 

0.82 

-0.52 

0.05 

Full fit with M g c 


Note. — The attenuation is the fraction of bolometric luminos¬ 
ity which is absorbed, averaged over all directions. 


1 Parameter in equation [7] 

2 The standard deviation of the scatter around the fit. 

3 This fit was performed with gas mass, instead of baryonic mass, 
driving the parameter <5 in equation [71 


simulations is 



r = T 0 


(1-0 
l ( L Y 

, where 

3 ( M* 

YY M b \ s 

0 

o 

v M 0 yr - 

- 1 / \ 10 11 M 0 J 


(7) 


and the Greek letters denote free parameters. 

The fits were done as simple x 2 minimizations, with 
one point for each simulation snapshot in the simula¬ 
tions listed in Table 0 As there were nine different Sbc 
merger simulations and only one of each of the simu¬ 
lations using the other galaxy models, the Sbc simula¬ 
tions were given only 1/9 the weight in the fits. This 
was done to avoid giving undue weight to the massive, 
bright, and gas-rich Sbc mergers. Because the time be¬ 
tween saved simulation snapshots was varied in order to 
capture short-lived stages, each simulation snapshot was 
also weighted in proportion to its “time of influence”, i.e. 
the time to preceding and following snapshots. Apart 
from this weighting, constant errors on the points were 
assumed. Table 0 shows the parameters resulting from 
fits of the bolometric attenuation to the simulations un¬ 
der different constraints. The quantities come in with 
different powers from those in equation El but, as noted 
earlier, this is not surprising. 

A plot of the actual attenuations against those pre¬ 
dicted from the full fit is shown in Figure 0 The fit 
describes the behaviour of the simulations well, with a 
la scatter of about 0.04 in the absorbed fraction. Since 
this fit describes the fraction of bolometric luminosity 
absorbed by dust, it can be used to predict the dust lu¬ 
minosity of these systems. It should be emphasized that 
the quantities used are those of the system, implying that 
when the merging galaxies are still distinct it is the ag¬ 
gregate luminosity, mass, and star-formation rate of the 
two which is used. 

Table 0 contains additional parameter sets besides the 
complete fit. These are appropriate if some of the quan¬ 
tities which go into equation 0 are unknown. With a 
modest increase in the scatter around the relation, one 
of the largely complimentary quantities L or M*. can be 
excluded. In these cases, the power of the quantity not 
excluded increases to assume the role of the excluded 
quantity. With significantly increased scatter, the at¬ 
tenuation can also be predicted from one quantity only. 
This is presumably because metallicity, luminosity, and 



Fig. 7. — Actual attenuation (fraction of bolometric luminosity 
absorbed by dust) in the simulations, compared to the attenuation 
predicted by the full fit of equation 0 (the first line of Table |5J. 
Each point is a simulation snapshot. 


mass are intrinsically correlated in galaxies. It is difficult 
to say to what degree these single-parameter correlations 
are affected by our small set of galaxy models. 

Finally, the last line of Table0contains the parameters 
obtained when the fit was done using the gas mass. The 
fit has about 50% larger scatter than when the baryonic 
mass was used, so it still provides a useful description of 
the data. However, the dependence on the quantities is 
strange, [3 is negative which means that more luminous 
galaxies should have smaller dust attenuation. It is as 
if luminosity has assumed part of the function of the 
baryonic mass. 

3.7. Fits At Specific Wavelengths 

The fits in the previous section are useful for predict¬ 
ing the infrared dust luminosity of the systems, but they 
cannot be used to predict the attenuation of radiation 
emerging from the system observed at a specific wave¬ 
length. Furthermore, in order to be able to make pre¬ 
dictions of the luminosity inferred by an observer, the 
variation in attenuation with line of sight must be con¬ 
sidered. Tabic 0 contains fits, also using equation 0 to 
the attenuation in each of the bandpasses included in the 
calculation. Here, each point in time of a simulation is 
associated with 11 different attenuations, one for each of 
the different viewing angles calculated. These fits have 
significantly larger scatter than the fit to the bolometric 
attenuation in Table 0 but this is largely due to the vari¬ 
ation of attenuation over different lines of sight. When 
the attenuation is averaged over all lines of sight, the 
fit is unchanged, but the scatter is much smaller, only 
marginally larger than for the fit to bolometric attenua¬ 
tion. An example of a fit, for the SDSS g band, is shown 
in Figure 0 

3.8. Observations of Dust Attenuation 

Several previous studies have attempted to estimate 
the dust attenuation in observed galaxies nature, and 
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TABLE 3 

Fits to the dust attenuation in specific filters. 


T0 a 

Q a 

/3 a 

7 a 

<5 a 

IT b 

&i C 

Filter 

7.67 

3.33 

-0.04 

0.29 

-0.74 

0.07 

0.04 

GALEX FUV 

5.03 

2.45 

0.03 

0.36 

-0.61 

0.07 

0.04 

GALEX NUV 

1.84 

1.52 

0.45 

0.28 

-0.70 

0.10 

0.05 

SDSS u 

1.27 

1.39 

0.55 

0.19 

-0.70 

0.11 

0.05 

SDSS g 

0.94 

1.36 

0.61 

0.16 

-0.71 

0.12 

0.05 

SDSS r 

0.78 

1.34 

0.65 

0.12 

-0.70 

0.11 

0.05 

SDSS i 

0.68 

1.37 

0.69 

0.10 

-0.72 

0.11 

0.05 

SDSS 2 

0.53 

1.43 

0.72 

0.09 

-0.74 

0.10 

0.05 

2MASS J 

0.42 

1.54 

0.76 

0.09 

-0.78 

0.09 

0.05 

2MASS H 

0.32 

1.62 

0.77 

0.10 

-0.80 

0.08 

0.05 

2MASS Ks 

0.19 

1.78 

0.74 

0.14 

-0.82 

0.06 

0.04 

Spitzer IRAC1 

0.14 

1.95 

0.71 

0.18 

-0.86 

0.06 

0.04 

Spitzer IRAC2 

23.9 


0.70 



0.12 

0.10 

1900A(WH96) d 

4.0 


0.49 



0.15 

0.11 

B (WH96) d 


Note. — Fits were performed using all parameters. The scatter 
in the far-ultraviolet bands is suppressed because most attenuations 
are close to 1. 


1 Parameter in equation |?| 

2 The standard deviation of the total scatter around the fit, including 
the variation with viewing angle. 

3 The standard deviation of the intrinsic scatter around the fit, i.e. 
excluding the variation with viewing angle. 

4 Fits done as a compar ison to the study by IWang &; Heckmanl 
J1996D . described in section m 



Fig. 8. — Actual attenuation in the simulations versus the atten¬ 
uation predicted by equation [71 but now for the luminosity in the 
Sloan Digitized Sky Survey g band. Unlike Figure [7] which shows 
absorbed energy averaged over all lines of sight, this figure shows 
the line-of-sight attenuation. Each point is a simulation snapshot 
from a certain line of sight. The slightly negative attenuations re¬ 
sult from preferential scattering out of the plane of the disk in the 
initial galaxies. This fit has significantly larger scatter than the 
fit to bolometric attenuation, a result mainly because of the large 
variation of attenuation with viewing angle, but provides a good 
description of the behavior averaged over all directions. 


o ur simulations can be com pared to these results. 

W’aii'j & Heckman (Il996t) examined the dust content 
in late-type galaxies based on their UV/FIR flux ratio 
and fit the derived optical depth to a power-law depen¬ 
dence on luminosity. To compare to their study, a fit 
of the ultraviolet attenuation using only a dependence 
on ultraviolet luminosity is also performed. This fit is 


shown in Figure |9j This fit was complicated by the fact 
that when using only luminosity instead of all parame¬ 
ters to predict the dust attenuation, there is a signifi¬ 
cant difference between the initial merging galaxies and 
the merger remnants. Especially in the smaller G2, Gl, 
and GO mergers, the merger remnants have significantly 
higher dust attenuation for their luminosity compared 
to the earlier stages. Furthermore, because the simu¬ 
lated galaxies spend a lot of time as merger remnants 
(limited by how long the simulation has been run), this 
stage has significant weight when fitting. These points 
can be seen in the upper left of Figure El As discussed 
below, our model likely overestimates the dust attenua¬ 
tion in the merger remnants. For this reason, and also 
because WH96 did not include early-type galaxies, these 
low-luminosity, high-attenuation points were excluded in 
this fit. (If these points are included, the resulting fit, 
also shown in Figure El does not provide a good de¬ 
scription for the low-luminosity pre-merger and merging 
galaxies.) 

After excluding these points, the fit yields (3 = 0.70, 
consistent with the result of WH96 who obtained (3 = 
0.5 ± 0.2. In terms of the normalization, our fit results 
in tq = 23.9, at a luminosity of 10 11 Lq. Rescaled to 
the WH96 luminosity zeropoint of 4.5 x 10 9 Lq, it cor¬ 
responds to a UV optical depth of 2.7, while their re¬ 
sult was 1.7 ± 0.6. The simulated galaxies seem to have 
larger optical depth for the same luminosity, possibly re¬ 
flecting the fact that they not really “normal” spirals. 
However, the WH96 galaxies were essentially limited to 
A1900 < 10 10 Lq, and in this luminosity range, the WH96 
fit is a good description also of the simulated galaxies. 
The intrinsic scatter in this fit is much larger than when 
all the parameters are used. 

WH96 also presented their results in terms of the 13- 
band luminosity and, correspondingly, the attenuation in 
the B band. To compare to these result, a fit to the simu- 
lations was d o ne in t he B band (using the prescription of 
iFukueita et alJlT996i to transform the SDSS magnitudes 
into a B magnitude). The same points were excluded 
for this fit, also shown in Figure El as for the ultraviolet 
fit. This results in [3 = 0.49, in excellent agreement with 
WH96. Their normalization is lower also in this case. 
Our result was r 0 = 4.0, which, rescaled to the WH96 
luminosity of 1.3 x 10 10 Lq, corresponds to an optical 
depth of 1.5. Their result was 0.8 ± 0.T_ 

A similar study was performed bv IViih et al.1 1 2003 1. 
who estimated the UV attenuation in a sample of Lyman- 
break galaxies using radiative-transfer models of clumpy 
dust shells. They basically obtained 

1 - ^ oc L ~°- 95±0 - 5 , (8) 

Jj 

at 1600 A. In optically thick situations, this would corre¬ 
spond to (3 = 0.95±0.5 in our formulation, also consistent 
with o ur resu lt at 1 900 A. 

iBuat et all 1 2004 1 compared the dust attenuation, in 
the GALEX NUV band, of two samples selected in the 
near-ultraviolet and far-infrared. While they saw an in¬ 
crease in dust attenuation with luminosity in both sam¬ 
ples, they did not attempt to determine the luminosity 
dependence. No attempt has been made to replicate their 
selection criteria in our simulated galaxies either, but the 
distribution of attenuation with luminosity in the sirnu- 





















Dust Attenuation in Simulations of Interacting Galaxies 


9 


O 

O 

CD 



CO 

-Q 

O 

_J 




.U _i_ ■ 3gs _,. .i _,. . 

10 8 10 9 10 10 1 0 11 1 o 12 


B/ L© 


Fig. 9.— Comparison with I Wang &: Heckmanl <119961) . On the left, luminosity at 1900 A versus attenuation at 1900 A for the simulations, 
on the right luminosity in the B band versus attenuation in the B band. The size of the symbols are proportional to the fitting weight of 
the points. The dashed (green in the Electronic Edition) lines show the fits marked “WH96” in Table l3l The points in the upper left re gion 
in the left plot, marked by the solid line, consist of merger remnants and were excluded from the fit for reasons explained in Section [3.71 
Because these points are not nicely delineated in the right plot, they are also marked as (red) crosses. (The dotted (red) lines show the fits 
obtained if these points are included in the fit.) The dot-dashed (blue) lines show the WH96 fits. In the ultraviolet, the WH96 fit agrees 
fairly well with the simulated galaxies for luminosities < 10 10 L©, the maximum UV luminosity of the WH96 galaxies, but underpredicts 
the attenuation at high luminosities. In the B band, the WH96 fit has a luminosity dependence similar to our fit, but is normalized lower. 
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Labs/ L© 


Fig. 10.— Luminosity absorbed by dust versus attenuation in 
the GALEX NUV band (in magnitudes) for the simul ated galax¬ 
ies. T his plot can be directly compared to Figure 3 of IBuat et alJ 
Eool^ . Their infrared-selected galaxy sample shows a very similar 
distribution, while their UV-selected sample has lower attenua¬ 
tions. Given that our simulated galaxies are bright star bursts, the 
infrared selection is likely more app ropriate, at least a t higher lu¬ 
minosities. For L abs < 10 10 L©, the lBuat et alJ 120041) IR-selected 
sample contains very few galaxies, probably a selection effect, but 
the ones present lie in the same region as the simulations. 


lations (shown in Figure HOI) is similar to that of their 
sample. 

3.9. Model Limitations 


While our suite of simulations is the most ambitious 
effort to self-consistently model dust extinction in star- 
bursts so far, our model has limitations. For one, dust is 
known to be clumpy or patchy on small scales, something 
which is not t aken into account in the simulations. Pre- 
vious studies llWitt fc Gordonlfibbfl l2f)0ft iGordon et alJ 
1112) have indicated that clumping has important im¬ 
plications for the effects of dust. While the large-scale 
structure of the dust distribution in the galaxies is de¬ 
termined by the hydrodynamic simulations, clumping on 
the scale of individual star-forming regions is far below 
our resolution and is not included. The most likely effect 
of clumping would be to open up “holes” with lower at¬ 
tenuation, so that our simulations would underestimate 
the frequency with which young stellar populations are 
visible. 

Another phenomenon which is treated poorly in our 
simulations is gas outflows. Starburst galaxies are al¬ 
most ubiq uitously observed to h ave large-scale gas out¬ 
flows (e.g. iHeckman et alJ 120001) . While our feedback 
prescription does provide energy input from supernovae, 
it does not seem to lead to significant outflows. Fur- 
therm ore, recent sim ul atio ns including active galactic nu¬ 
clei ( ISpringel et al.l[2005btl show that AGN feedback can 
drive a powerful outflow, clearing out gas from the cen¬ 
tral regions of the galaxies. This phenomenon is also not 
included in our present simulations. In the context of 
dust attenuation, this lack of outflows has two implica¬ 
tions. First, outflows, like dust patchiness, would tend to 
open up holes in the dust distribution. Second, the lack 
of outflows also means that the metals produced by the 
starburst remain in the starburst region. Since the op¬ 
tical depth is determined using a constant dust-to-metal 
ratio, this will translate to a larger opacity within the 
star-forming regions than would be expected if a signifi- 
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cant fraction of the metals were ejected. 

Both of these limitations make it probable that our 
simulations overestimate the amount of dust attenua¬ 
tion to some degree, especially in the late stages after 
the galactic nuclei have merged. Less clear, however, is 
whether the correlations found here would be ruined. For 
example, while an AGN-driven outflow will clear out gas 
and dust from the galaxy and hence decrease the dust at¬ 
tenuation, it also truncates star formation. Thus, while 
the evolution of the (stellar) luminosity of the starburst 
would be altered by this phenomenon, the correlation 
between dust attenuation and the luminosity, metallic- 
ity, etc., of the starburst would not necessarily be al¬ 
tered. Work is underway to improve our model and in¬ 
clude these effects. 

4. DISCUSSION 

One of the remarkable aspects of Figure El is that the 
low-luminosity stages of the simulations, which comprise 
both the initial, separate spiral galaxies and the final 
merger remnant, overlap. A priori, there seems to be 
no reason to expect that the attenuation should remain 
constant if two galaxies are merged into one. However, 
from the fit results in Table 0 it is evident that the pow¬ 
ers of the extensive quantities L , M*, and Mb approxi¬ 
mately add to zero. This results in a dust attenuation 
which is almost insensitive to simple size scaling of the 
galaxies; rather than being determined by luminosity or 
star-formation rate alone, the dust attenuation seems to 
be governed by “specific luminosity” and “specific star- 
formation rate”, i.e., L/Mb and This is fortu¬ 

nate for our analysis, since it means that not treating the 
two galaxies in the initial stages of the merger separately 
does not bias the results. This is only true in the spe¬ 
cific case being treated here, where the two galaxies are 
identical. If the two galaxies were not identical, as when 
simulating minor mergers or a merger between a spiral 
and an elliptical galaxy, the dust attenuation would have 
to be determined separately for the two components re¬ 
gardless of whether a, (3, and 7 sum to zero. 

It should, however, be pointed out that a + (3 + 7 « 0 
is a poor approximation for the fits to the attenuation 
in the GALEX bands in Table [3] There is thus the pos¬ 
sibility that the fits in these bands have been biased by 
the use of system, rather than individual galaxy, quan¬ 
tities. Indeed, there is a significant discrepancy between 
the fits and the initial stages of the simulations in these 
bands, such that the fits overestimate the amount of dust 
attenuation. 

Looking at the parameter sets in Table 0 clear trends 
with wavelength can be seen. Going from ultraviolet 
to near-infrared wavelengths, the dependence on star- 
formation rate decreases while the dependence on lumi¬ 
nosity increases. This is not unexpected; the luminosity 
in the ultraviolet is dominated by massive, short-lived 
stars, and hence correlates well with star-formation rate. 
At longer wavelengths, contributions to the luminosity 
come from stars of a wide range in age and is better rep¬ 
resented by the bolometric luminosity of the galaxy. This 
trend with wavelength thus contains information about 
the stars whose radiation is being absorbed. What is 
more surprising is that the trend is reversed at wave¬ 
lengths longer than 2 /im. Naively, one would expect 
that the dust attenuation at progressively longer wave¬ 


lengths always would be more dominated by older stars, 
but this does not seem to be the case. This effect proba¬ 
bly originates in the fact that around the age of 10 Myr, 
a stellar population is very bright in the near-infrared 
due to the presence of red supergiants. This means that 
young stellar populations make a larger contribution to 
the emission at several pm than what is expected from 
their main-sequence temperatures. 

It was earlier noted that, as shown in Figure 0 the 
UV/visual luminosity not absorbed by dust essentially is 
independent of the intrinsic bolometric luminosity. This 
notion is confirmed by the fit to the bolometric attenu¬ 
ation independent of the star-formation rate in Table 0 
That fit yielded (3 = 1.10, close to (3 = 1 for which the 
increase in luminosity is exactly compensated by the in¬ 
crease in attenuation given by equation 0 (in the optically 
thick limit). 

Is the fact that the dust attenuation in the simula¬ 
tions is well described by equation 0 indicating something 
about the relative geometry of dust and stars in the simu¬ 
lations? As already mentioned, gas density is the driving 
factor behind both dust optical depth and star-formation 
rate. Furthermore, in the simulations, dust and stars are 
assumed to be uniformly mixed within individual grid 
cells. It is thus not unreasonable to expect that a uni¬ 
form mixture of dust and stars should fit the simulations 
reasonably well. 

The fact that the dust attenuation can be predicted 
so well by a simple formula also implies that there are 
no hidden parameters determi ning whether a luminou s 
galaxy will be a ULIRG or not. iBekki fc Shiovai (l2000bl) . 
from analyzing a prograde-prograde and a retrograde- 
retrograde major merger at the time of highest star- 
formation rate, drew the conclusion that retrograde 
mergers should have stronger internal dust attenuation 
than prograde ones. This led them to conjecture that in¬ 
teracting galaxies without long tidal tails should be more 
prevalent among ULIRGs. This result is not confirmed 
by our much more extensive analysis. While there is a 
tendency for retrograde mergers to induce slightly more 
intense starbursts in our simulations, and hence be more 
obscured according to our fitting formula, they follow 
the same relation for dust attenuation as mergers of any 
other geometry. At least within the parameter space cov¬ 
ered by these simulations, there should be no such thing 
as a “naked” vigorously star-forming, 10 12 L e system. 
While the attenuation along some lines of sight might be 
smaller, the vast majority of the bolometric luminosity 
should always be emerging as infrared dust emiss i on. 

Com paring the s tudies of IWang fc Heckman! ( 1996 1 
a.nd lViih et al . ( 1200 , 1 I f to these simulation results, it is en¬ 
couraging to note that the dependence of optical depth 
on luminosity is similar. Unfortunately, they did not 
have metallicity or mass information, so they were un¬ 
able to perform a multi-dimensional fit like the one in 
this study. It would be interesting to know whether the 
dependence on these parameters would also be similar to 
what is found here. 

The fitting formulae presented here are in the “theoret¬ 
ical plane”, i.e., they depend on quantities such as total 
baryonic mass and bolometric luminosity which are gen¬ 
erally known in theoretical models but hard to determine 
from observations. As such, they are of limited use for in¬ 
terpreting observations. It has already been emphasized 
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that dust effects cancel the effects of increased luminos¬ 
ity to a remarkable degree, so that the simulations show 
virtually no correlation between the apparent luminos¬ 
ity and the dust attenuation, either bolometric or in the 
ultraviolet. Unfortunately, this implies that drawing con¬ 
clusions about these highly dust-extinguished systems 
from their apparent luminosity is very difficult, unless 
infrared data are available. In particular, “correcting” 
apparent luminosities for dust using a luminosity-based 
pr escri p tion is no t likely to work well, a fact also noted 
bv lHookins et all 1200 ll) . 

The relations presented here should be suitable for in¬ 
clusion in theoretical models for galaxy formation, such 
as semi-analytic models (SAMs) includin g merger-driven 
starbursts fe.g.. iSomerville et al .1 120011) . Unlike other 
theoretical models fo r dust effects (e.g.- lSilva et al.lfl998t 
iCharlot fc Fallll200(Tl our results give information about 
the magnitude of the dust attenuation and its depen¬ 
dence on the properties of the galaxy in a self-consistent 
way. Current approaches used to incorporate dust in 
SAMs include relying on empirical results like those of 
WH96 or even simpler approximations such as a uni¬ 
form slab with optical depth proportional to gas column 
density times metallicity (essentially the toy model pre¬ 
sented h ere). Even in more so phisticated models, such as 
those of lGranato et al.1 (I2000I) . the optical depth is typi¬ 
cally determined in the same simple way. These models 
also contain numerous adjustable parameters whose val¬ 
ues are not given by the SAM and thus must be fixed at 
constant values. Given that our fits appear to work well 
across a wide range of galaxy properties, including these 
results in SAMs should lead to a more realistic estimation 
of the effects of dust in cosmological scenarios. 

Finally, nothing has been said here about the shape of 
the dust attenuation cu rve, e.g., if the simula ted galaxies 
obey the “Calzetti law” llCalzetti et a,l.l1994lh This infor¬ 
mation is contained in our fits, and a future paper will ex¬ 
plore this in detail. For now, we simply remark that our 
simulations are inconsistent with a simple “screen-like” 
attenuation curve like the Calzetti law, and do not even 
follow any single red dening law resulting from more com¬ 
plicat ed geometries flWitt fc Gordonll996t iGordon et all 
EH). This is likely the result of a stochastic superpo¬ 
sition of many star-forming regions with various optical 
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depths and ages. 


5. SUMMARY 

We have presented results from radiative-transfer cal¬ 
culations in a comprehensive suite of galaxy merger sim¬ 
ulations. The results from these simulations consist of 
images at many different wavelengths, as well as spec¬ 
tral energy distributions, and seem to agree well with 
observed properties of starburst galaxies. 

It was discovered that the dust attenuation (defined 
as the fraction of luminosity which is absorbed by dust) 
in the simulations can be predicted from the bolomet¬ 
ric luminosity, star-formation rate, baryonic mass, and 
average gas metallicity of the system through a simple, 
physically motivated formula (eq. 0. Averaged over all 
directions, the attenuation of the bolometric luminosity 
can be predicted with a scatter of 4%. The attenuation 
along a specific line of sight can be predicted with a scat¬ 
ter of 6 - 12%, depending on wavelength. The increased 
scatter is largely a result of the variation of attenuation 
with viewing angle. These relations are valid for simula¬ 
tions with a range of two orders of magnitude in mass, 
with metallicities from 0.3 to 1.1 Zq, and gas fractions 
from 20 to 60%. The relations also seem to be valid for 
both isolated and interacting galaxies. 

Our result s are consistent w i th stu dies o f observed 
galaxi es by iWang & Heckmanl 1199611 and iViih et all 
1 2003 ). but our inclusion of additional independent vari¬ 
ables significantly lowers the scatter around the relations, 
and we present results valid for wavelengths from the far- 
ultraviolet to the near-infrared. 
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